Calculation of hysteretic force between a superconductor and a magnet 
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The magnetic levitation forces of a superconducting disk (SC) levitated by a cylindrical permanent 
magnet (PM) have been calculated from first principles for superconductors with finite thickness. 
The current j{p,z) and field B{p,z) profiles in the SC in the non-uniform magnetic field generated 
by the PM have been derived. The levitation force is observed to depend non-linearly on the critical 
current density jc and on the thickness of the SC. The flux creep is described by a current-voltage 
law E{j) — Ec{j / jc)" , from which we show that the levitation force depends on the speed at which 
the PM approaches or recedes from the SC, which accounts for the experimentally observed force 
creep phenomenon. The stiffness of the system has been derived by calculating the minor force 
loops. The numerical results reproduce many of the features exhibited by experiments. 

PACS numbers: 74.60.-w, 74.25.Ha, 74.25.Ld 



I. INTRODUCTION 

It has been well known that a permanent magnet (PM) 
can be stably levitated above a high temperature su- 
perconductor (HTS) cooled by liquid nitrogen, which 
has become the well known symbol for HTS technol- 
ogy. This fascinating magnetic levitation results from the 
interaction of the induced current inside the supercon- 
ductor with the inhomogeneous magnetic field generated 
by the PM. Because of its possible industrial appli|ea|- 
tions, such as noncontacted superconducting beaxinjss- 
gravimetersEl, flywheel energy storagej-systemsoLiiJiD'Effi , 
magnetic levitation transport systemsE^, and motorsHil, 
the magnetic levitation between a PM and a HTS has 
been the subject of intensive studies for the last decade. 

The most common feature of the magnetic levita- 
tion is the hysteretic behavior of the vertical force 
versus the distance z between the PM and the HTS 
when the PM is descending to and then ascending from 
the HTS. When a PM approaches a zero- field-cooled 
HTS, the levitation force increases monotonically from 
zero, as the PM is moving from the HTS, the levita- 
tion force decreases sharply to a negative peak at some 
distance, indicating attractive force between the HTS 
and the PM, then declines to zero again at larger dis- 
tance. Detailed experiments have been performed on 
melt-textured-grown (MTG) YBCO and RBCO (R de- 
notes a rare earth element such as Nd, Sm, Gd, Eu, 
Dy, Ho, Er, Tm, Y.h, .T;ij_and La), thin films, as well 
as granular samplesyBOldEl. For MTG YBCO and 
RBCO samples, the force curve is usually asymmetrical, 
i.e., the absolute value of the attractive force is smaller 
than the maximum repulsive force, |i^a,max| < ^r,max- 
For thin films, tiie force curve is almost symmetrical 
|-Fa,max| ~ Jr.maxli^, whilc for [granular samples, the at- 
tractive force is hardly presentll3. 

Various parameters are very important in determining 



the levitation force. The most important are the pinning 
strength (critical current density) of the HTS and the in- 
duced shielding current loops inside the HTS. Higher jc 
and larger loops are very important to achieve a high lev- 
itation force. It has been suggested that the levitation 
force increases linearlsi-|With jc and more than linearly 
with the current loopsEJ. According to this criterion, the 
current material of choice for superconducting levitation 
is MTG YBCO and RBCO samples, because high quality 
MTG samples with high critical current density and sin- 
gle domain diameters up ttiilO cm are currently available 
by means of melt procesalj. Another reason is its high 
irreversibility line at liquid nitrogen temperature, which 
is the normal operating temperature for the magnetic 
levitation experiments. Although higher critical current 
density can be achieved in thin films, the levitation force 
is limited by the thickness of the films. 

According to the critical state modelElH, the levita- 
tion force is independent of the speed at, which the PM 
approaches and recedes from the HTSc^. However, as 
the thermally activated flux motion is prominent in HTS, 
resultingpiii the relaxation of the magnetization (current 
density) ,c3 the levitation force which depends on the cur- 
rent density of the HTS is expected to decrease with time 
(force creep). Experimentally the levitation force is ob- 
served to be approximately logarithmic in time and can 
be well correlated with the thermally activated flux mo- 
tion in HTSB. 

The thickness of the HTS also drastically influences 
the levitation force. Because the critical current density 
is limited, for thin HTS samples, the levitation force in- 
creases linearly with the thickness of the HTS. However, 
beyond a certain thickness, the levitation force is inde- 
pendent of the thickncsstJ. Other characteristics of the 
HTS, such as the anisotropy and the grain orientation 
inside the MTG samples lia,Y6,been studied and shown to 
affect the levitation forceEJta. The stiffness of the PM- 
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HTS system, which represents the spring constant asso- 
ciated with vibrational motion of the levitation system 
has been intensively studied. And Experimental results 
for vertical stiffnesSj-ilateral stiffness and cross stiffness 
have been prescntcdolla. 

Besides these, the geometry and condition of the PM 
also affect the levitation force. The size and shape of the 
PM, the homogeneity and the temperature dependence 
of the PM magnetisation have been shown to influence 
the levitation forcet3. 

Although the magnetic levitation between a PM and a 
HTS can be easily demonstrated and detailed experimen- 
tal results for the levitation force have been presented, 
the theoretical models for this magnetic levitation sys- 
tem have not been fully developed yet so far. The first 
reason is that a model has to consider a finite thick- 
ness superconducting disk under a perpendicular mag- 
netic field. In the case of an infinite long cylinder under 
a parallel applied magnetic field, the magnetic field and 
current density profiles can be easily obtained by means 
of the critical state model. However, for a superconduct- 
ing disk under a perpendicularly applied magnetic field, 
the extreme demagnetization effects make this case qual- 
itatively different from the parallel field case. The second 
reason is that the model has to consider the response of a 
HTS immersed in a non- uniform magnetic field generated 
by the PM. 

In order to avoid the above difficulties, all the models 
for the magnetic levitation force presented so far have 
taken different assumptions based on different physical 
considerations. . Mos t of the models are based on the crit- 
ical state modelli3E3, which has long been used to account 
for the irreversiblft^paperties of type-II superconductors. 
The image modelcBO'Elreats the superconductor as a 
pure diamagnet and the PM as a set of magnetic dipoles, 
if the hysteretic behavior is not important. The system is 
described by two magnetic dipoles, representing the PM 
and its mirror image, oppositely magnetized and located 
the same distance below and above the superconductor 
surfaces. The advantage of this model is that analytical 
expressions can be obtained. However it cannot be used 
to account for the dynamic stiffness. The extension of 
the image model to include another 'frozen' image of the 
PM has beeniHtroduced to explain the dynamic stiffness 
of the systemcj. 

In a recent article by Navau and SanchezES, the authors 
reviewed the model in the critical state. Early model con- 
sidered only extreme limits such as complete flux exclu- 
sion or complete fiux penetration, thus described the be- 
havior of type-I superconductors or. tyne-II superconduc- 
tors with very high critical currentoO. This modeL 
later extended to describe the fiux penetratioi£3Sl 
However, these models assumed a superconducting sam- 
ple small enough to consider the magnetic-field-gradient 
constant along it. At the same time, the demagnetiza- 
tion effect due to the finite dimensions of the supercon- 
ductor has been neglected. Navau and SanchezE2l have 
accounted for the demagnetization by introducing a de- 



magnetization factor. However, it should be pointed out 
that the use of demagnetization factor for superconduc- 
tors in the mixed state in which current is distributed 
inside the bulk sample, is qualitatively invalidHtj. Al- 
though the authors considered a non-uniform field gradi- 
ent along the sample, the radial magnetic field generated 
by the PM has been neglected. While in real experi- 
ments, the PM is usually smaller than the HTS, in this 
case the assumption of neglecting the radial field may not 
be applicable. 

The models mentioned above can be used to account 
for some of the levitation force features exhibited by ex- 
periments, however, they are obviously oversimplified. 
First, the current and field distribution in the HTS under 
a non-uniform magnetic field have been neglected. For 
practical applications using the magnetic levitation, the 
precise control of the interacting force between a super- 
conducting target and the applied non-uniform magnetic 
field is of critical importance. It is hence necessary to 
thoroughly understand the hysteretic force characteris- 
tics of a superconductor immersed in such fields. From 
an analytical point of view, a precise knowledge of cur- 
rent and field profiles is essential to the force calculation 
of the magnetic levitation force between the PM and the 
HTS. Second, the models are all based on the critical 
state model, which neglects flux creep effects. However, 
at the nitrogen temperature, which is the normal operat- 
ing temperature for the levitation force experiments, the 
fiux creep may play an important role on the levitation 
force. Numerical analysis such as finite element method 
(FEM) has been applied to calculate the magnetic levi- 
tation force, it is capable of deriving the current and field 
profiles inside the HTS, but it neglects the fiux creep as 
well. 

To better understand the design of the magnetic levita- 
tion system between the PM and the HTS, a better model 
is required. In this paper, taking all the above points 
into account, we develop a model to describe a more re- 
alistic experimental configuration and to account for the 
experimentally observed features of the magnetic levita- 
tion force. The demagnetization effect is considered by 
calculating the current and magnetization of finite thick- 
ness superconductors, and current and field distributions 
in the superconductor are derived. The flux creep ef- 
fect is taken into account by using a current-voltage law 

Eij) = E,ij/j,y\ 

The paper is structured as follows. In Sec. II we will 
discuss the theoretical consideration for both the perma- 
nent magnet and the superconductor in a non-uniform 
magnetic field, including assumptions made and numer- 
ical considerations. The results of the calculations, in- 
cluding the effects of the HTS geometry and character- 
istics on the levitation force, the force creep effects, the 
minor loops and the stiffnesses will be presented and com- 
pared with experiments in the literature in Sec. Ill, and 



Sec. IV are the concluding remarks. 
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II. MODELING 
A. The Permanent magnet 

We consider a superconducting disk with radius a and 
thickness 26, levitated over a co-axial cylindrical perma- 
nent magnet with radius i?pM and thickness ipM- The 
top surface center of the PM is taken as the origin of the 
cylindrical coordinate system (p, 0, z). Because of the ax- 



ial symmetry of the system, only the cross section of the 
system is considered, with the axis direction Z chosen as 
the symmetry axis of the SC and the PM, and p (radial 
direction) parallel to the surfaces of the SC and the PM. 
For this configuration, the vector potential of the PM 
has only one component along the (j) direction and can 
be derived by integrating the vector potential of a circu- 
lar current loop with radius i?pM along the thickness tpM 
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where Brom is the remanent induction of the PM. The 
radial field Bp — —dA^/dz can then be written as. 
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B. The superconducting disk in a non-uniform 
magnetic field 



The response of superconducting disks and strips in 
an uniform applicc^jpagnetic field has been extensively 
studied by BrandtEZlES. We now consider the response 
of a superconducting disk in the non-uniform magnetic 
field generated by the PM [Eqs. (|) and §)]. The cal- 
culation is for the zero-field-cooled (ZFC) process, in 
which the superconducting disk is cooled below its crit- 
ical temperature Tc in zero applied field and then the 
non-uniform field is applied to the disk. Screening cur- 
rent is generated inside the superconducting disk. The 
central idea is to find the equation of motion for the 



current density. Because of the axial symmetry, the 
current density J and the vector potential Aj gener- 
ated by the current have only one component along the 
(j) direction. The total vector potential of the system 
is then A = Aj + A^, and the total magnetic field 
_B = VxA = Vx Aj. We assume here the material 
law B — /ioH, which is a good approximation when the 
flux density B and the critical sheet current JJ) are larger 
than the lower critical field Bd everywhere inside the su- 
perconducting disk. This requirement is often satisfied in 
the magnetic levitation measurement, normally operated 
at a relatively high temperature of 77 K. According to 
the Maxwell equation with gauge V ■ Aj — 0, we have 
J = VxB = VxVxA = VxVxAj = -V^Aj. The 
solution of this Laplace equation in cylindrical geometry 
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can be written as 



Aj{p,z)^~ii^ dp' dz'Q{r,r')J{r'), (4) 

Jq J-b 

with r = (p, z) and = (p', z'). The integral kernel 
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if and E are the complete elliptic integrals of the first 
and second kind, respectively. 

Eq. (H) is obtained by integrating the 3D Green func- 
tion of the Laplace equation, l/47r|r3 — jp4| with r-^ = 
{x,y,z), over the angle (j) = arctan(y/a;)e3. The total 
vector potential can then be written as 

pa pb 

A{p,z) = ~pQ dp' dz'Q{r,r')J{r') + A^{p,z) 

Jo J~b 

(7) 

To obtain the desired equation of motion for the 
current density J{p,z,t), we express the induction law 
V X E = -B = V X A in the form E = —A. Combining 
this expression with Eq. (Q), we have 



by Eqs. (g) and (g), but one may also use other type of 
magnetic field, such as the field generated by a circular 
current loop or by dipoles. Besides these, the equation 
also depends on the material law E — E{J) of the super- 
conducting sample. Obviously any sufficiently sharply 
bent E{J) law may be used; in the following we choose 
the rather general material law E{J) = Ec{J/ Jc)"', which 
yields the limits of the critical state model for n —>■ 00, 
the flux flow model for n = 1, and the flux creep model 
for 1 < n < 00. Eq. is easily time integrated by 
starting with J{p,z,t = 0) = at time t = and then 
putting, J(p, z,t ^ t + dt) = J(p, z, t) + J{p, z, t)dt. The 
vector potential can then be derived from Eq. (|^) and the 
magnetic field is i? = V x A. For all the calculations in 
this paper, we use a reduced unit of i?c = a = po = 1- 
And the critical current is assumed to be independent of 
the magnetic field. 



C. The levitation force and the hysteresis loop 

The PM approaches and recedes from the HTS as 
z = Zf) — zo sin(cjt) + Zoo, where zo -I- zoo is the initial 
distance and zoo is the minimum distance between the 
PM and the HTS, lu represents the speed at which the 
PM approaches and recedes from the HTS. Experimen- 
tally, uncertainty will be caused when the PM touches the 
HTS, and therefore the limit zoo = should be avoided. 
In this calculation we choose zoo/a = 0.1 as the minimum 
distance between the PM and the superconducting disk. 

As the current density and the magnetic field inside 
the superconducting disk have been derived, the vertical 
levitation force along the z-axis can be readily obtained 
as 



E{p,z)^po dp' dz'Q{r,r')J{r')-A4p,z) (8) 

Jo J-b 

The equation of motion for the current density can be 
derived by inverting Eq. (||) as 

z) = — r dp' f dz'Q-\r, r')[E{r, r') + A^{p' , z')], 

(9) 
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where Q ^ is the reciprocal kernel defined by 
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It can be seen from Eq. (^) that the equation of mo- 
tion for the current density contains the reciprocal kernel, 
which depends only on the geometry of the superconduct- 
ing sample. Here we consider a disk, but the kernel Q can 
be changed to calculate other sample geonaetries, such as 
strips, and even arbitrary sample shapesE3. The equa- 
tion also depends on the applied magnetic field via A^. 
Here we consider the field of a permanent magnet given 



dp / dzJ{p, z)Bp{p, z) 



(11) 



As the applied magnetic field is non-uniform, in order 
to show the magnetization hysteresis loop, we choose the 
axis field -Baxis ^ Bz{p — 0, z) at the bottom surface cen- 
ter of the superconducting disk to represent the strength 
of the applied non-uniform field, with 



and the magnetic moment 



m — 2tt dp dzp J{p, z) 
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(12) 
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D. The minor force loop and the stiffness 

Because of the hysteretic behavior of the superconduct- 
ing disk, when the PM is moved away and back again by 
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z/a=0.48 



FIG. 1: Magnetic field lines when a PM approaches a su- 
perconducting disk with b/a = 0.25, Jc = 1, and o" = 20 at 
distances z/a — 1.79,1.19,0.92,0.48. Dashed lines are the 
contour lines of the current density inside the disk. The PM 
approaches the disk as z = zq — zq sm{ujt) with a frequency of 

LO = 0.1. 

a small amount 5z at a distance z, the Fz{z) curve follows 
a minor force loop rather than the major force loop. For 
small Sz, the minor force loop is reversible, but beyond 
a certain Sz, hysteretic behavior in the minor force loop 
will be observed. The stiffness at different distance z can 
then be derived as — —dF^/dz. 

III. RESULTS AND DISCUSSIONS 
A. Current and field profiles in the HTS 

The magnetic flux lines when the PM is approaching 
the superconducting disk with side ratio b/a — 0.25 and 
cr = 20, Jc = I are shown in Fig. |l| as solid lines at 
distances z/a = 1.79, 1.19, 0.92 and 0.48. The dashed 
lines are the contour lines of the current density inside 
the disk. Fig. |^ shows the penetration of the non-uniform 
magnetic field generated by the PM into the supercon- 



ducting disk. Some features different from those of the 
superconducting disk in uniform magnetic fieldEZI can be 
clearly seen. Because the magnetic field is stronger at 
the bottom of the disk (see z/a — 0.48 of Fig. |l|), the 
penetration starts from the bottom of the disk, while the 
top surface is not penetrated. The resulting magnetic 
field fronts inside the disk form an onion-shape rather 
than the symmetrical lens-shape observed for homoge- 
neous applied field. 

Fig. 1^ shows the profiles of the current density J{x, y) 
corresponding to Fig. |l|. The current density first satu- 
rates at the edges of the bottom surface, while the cur- 
rent density at the edges of the top surface is smaller 
{z/a — 1.79). As the PM is moving closer, the satura- 
tion spreads both into the top surface and into the middle 
of the disk {z/a = 1.19 and z/a — 0.92), until it is satu- 
rated everywhere inside the disk {z/a — 0.48). Because 
we choose a = 20, the maximum current density is much 
smaller than the critical current density Jc = 1. Increas- 
ing cr will result in a larger current density close to Jc- 

The magnetic field lines when the PM is moving away 
from the superconducting disk with side ratio b/a — 0.25 
and a — 20, Jc — 1 are shown in Fig. || as solid lines at 
distances z/a = 0.12, 0.20, 0.32 and 0.48. The dashed 
lines are the contour lines of the current density inside 
the disk. And Fig. ^ shows the corresponding profiles 
of the current density J{x,y). As the magnetic field is 
decreased, the current density is reversed inside the disk. 
What is interesting is that the reversion begins at the 
bottom surface with x/a = ±0.5, where the magnetic 
field is strongest. As the PM is moving away further, 
the reversion spreads onto the top surface and middle 
of the disk (Fig. |, z/a = 0.20 and z/a = 0.32), until 
the current is completely reversed {z/a > 0.48). Note 
that for a superconducting disk immersed in a uniform 
magnetic field, the reversion usually begins at the edge 
of the disk and spreads into the center of the sample. As 
the levitation force is determined by the current density 
J{x, y) and the magnetic field shown in Fig. |l| to Figjj, 
the features shown in these figures will be reflected in the 
levitation force as will be discussed in the following. 

The magnetic field lines and the profiles of the current 
density J(x, y) depend on the side ration h/ a, a and crit- 
ical current density Jc of the superconducting disk, as 
well as the dimensions of the PM. Detailed results will 
be presented elsewhere. 



B. Geometry of tlie HTS on 

As we use reduced units of -Be = a = Mo = 1 in this 
calculation, the effect of the HTS geometry on the levi- 
tation is demonstrated by calculating the levitation force 
for different thicknesses 26 of the superconducting disk. 
The results are shown in Fig. ||a, where ^7(2^) curves at 
a — 2, uj = 0.1 and Jc — ^ for different thicknesses 26 
are plotted. The F^, curves show typical hysteretic be- 
havior. For larger side ratio h/a, the hysteretic loop is 
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FIG. 2: Current profiles J{x, y) 

asymmetrical, with the maximum repulsive force larger 
than the maximum attractive force. However, as the 
side ratio is decreased, the hysteretic force loop becomes 
more and more symmetrical. For h/a = 0.05, the curve 
is completely symmetrical. Experimentally, symmetrical 
Fz -jiersus z curves have been observed in YBCO thin 
filmtj, while in MTG bulk samples, versus z curves 
are usually asymmetrical. The reason for this is easy to 
understand: For a thin sample, the magnetic moment is 
saturated when the PM is close to the sample. When the 
PM is moving away from the sample, only a little decrease 
in the applied field saturates the magnetic moment in the 
reverse direction, results in a symmetrical magnetization 
hysteresis loop (see Fig. where the magnetization hys- 
teresis loops corresponding to Fig. |^ are plotted), and 
therefore a symmetrical force loop. While for a bulk sam- 
ple the magnetic moment saturates when the PM is close 
to the sample, it is never saturated in the reverse di- 
rection when the PM is moving far away from the disk 




for the same disk as in Fig. |l[ 

(see the magnetization hysteresis loops for h/a > 0.8 in 
Fig. ^), resulting in an asymmetric hysteresis loop of 
the magnetization and therefore an asymmetric loop of 
the force versus distance. 

Another interesting feature shown in Fig. ^ is that a 
maximum is found in the repulsive force when the PM 
is at some distance away from the HTS, rather than at 
the minimum distance. Exp epiment ally this maximum 
has been observed in thin filmsEJ, while it was hardly ob- 
served in MTG samples. Riise et al. accredited this max- 
imum to the dimensions of the PM, which tends to vanish 
with larger tpM and smaller RpM- And the disappearance 
of the maximum in MTG sample was explained noting 
that bulk sMnples are less sensitive to a non-uniform field 
than filmaisL Another explanation comes from Sanchez 
and Navaiuj, who claimed that this maximum is a result 
of the minimum in the derivative of the field produced 
by the PM, based on a constant-field-gradient modelE^. 
They concluded that the maximum depends crucially on 
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z/a=0.48 



FIG. 3: Magnetic field lines when a PM is moving away from 
a superconducting disk with hja — 0.25, Jc = 1, and cr = 20 
at distances z/a = 0.12, 0.20, 0.32, 0.48. Dashed lines are the 
contour lines of the current density inside the disk. The PM 
recedes from the disk as z = zo — zo sin(ix)i) with a frequency 
of tj = 0.1. 



the side ratio b/a, the larger the side ratio b/a, the closer 
the maximum shifts to zero PM-HTS separation. When 
b/a is sufBciently large, the maximum may not exist. 
However, as can be seen from Fig. the maximum is 
independent of the side ratio b/a. The peaks are ob- 
served to be at the same distance for all side ratios and 
even for 6 S> a not shown here. In a separate calculation, 
we have found that this maximum is independent of the 
dimensions of the PM. We conclude that the peak arises 
from the intrinsic properties of the HTS. As will be seen 
from figures below, the maximum depends on the creep 
exponent a. The larger cr, the closer the maximum shifts 
to zero PM-HTS separation. When a is sufficiently large, 
the maximum may not be observed any more. 

Because a — n — 1, and n can be related to thCj^jp.- 
ning potential of the sample as n = UoiT, B) / kBTEIic3, 
smaller a means lower pinning potential or higher tem- 
perature. Both MTG YBCO and YBCO thin film have a 
relatively high pinning potential, however experimentally 



when a PM approaches a thin film, it may increase the 
local temperature on the film, resulting in a lower a, and 
therefore in the maximum. On the other hand, although 
the PM may increase the temperature of the surface of 
the bulk sample, the currents flow in a much larger vol- 
ume and the levitation force is determined by the bulk 
properties, therefore the maximum may not be observed 
in MTG bulk samples. 

In the inset of Fig. |^a, we show the maximum repulsive 
force as a function of the side ratio b/a. It can be seen 
from the figure that for small side ratio b/a, the maxi- 
mum repulsive force increases linearly with b/a, but sat- 
urates as b/a is further increased. Technically, a super- 
conducting disk with diameter 2a approximately equal 
to the thickness may be optimum for magnetic levita- 
tion, further increase of the thickness will not enhance 
the levitation force. This calculated, result is consistent 
with the experimental observationfj. 



C. Characteristics of the HTS on Fz 

Characteristics of the HTS in this calculation are rep- 
resented by two parameters. One is the creep expo- 
nent a related to the pinning potential by n = cr -I- 1 = 
Uo{T, H)/kBT; another one is the critical current density 
Jc representing the pinning strength. The effect of a on 
the magnetic levitation force is shown in Fig. where 
the vertical magnetic levitation force versus the dis- 
tance z at b/a — 1.4, u; ~ 0.1 and Jc ~ 1 for different 
pinning potentials of the superconducting disk are plot- 
ted. It can be seen from Fig. ^ that as a is increased, the 
hysteretic force loop increases and the maximum repul- 
sive force shifts to smaller PM-HTS separation. When 
a is sufficiently large, the peak is not observed, as dis- 
cussed in the above section (the maximum at a — 100 
corresponds to minimum HM-HTS separation). On the 
contrary, the maximum in the attractive force shifts to 
larger PM-HTS separation with increasing a. The max- 
imum attractive force increases with cr, however, when a 
is larger than 5, it decreases again and saturates at higher 
a. As opposed to this, the maximum repulsive force in- 
creases monotonically with a and reaches a saturation 
value at large a, as shown in the inset of Fig. ||a. 

Fig. ^ shows the corresponding magnetization hys- 
teresis loops. For (7 = 1, the pinning potential is very 
small, because of the relaxation effects, the applied mag- 
netic field can penetrate deeper into the sample, similar 
to the case of small side ratio b/a shown in Fig. ^. The 
magnetization hysteresis loop and therefore the force loop 
are symmetric. Increasing the pinning potential (larger 
a) results in larger hysteresis loops of the magnetization. 
However, further increase of a leads to a reversible re- 
sponse and to a smaller magnetization hysteresis loop. 
Interestingly, the peak in the maximum attractive force 
corresponds to the maximum hysteresis of the magneti- 
zation (here cr = 5) as can be seen from Fig. ^ which is 
also the case in Fig. ^. 
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FIG. 4: Current profiles J{x, y) 

The sharp decrease of the levitation force from repul- 
sive to attractive results from the non-uniform magnetic 
field generated by the PM and from the relaxation of the 
current in the HTS. According to Eq. (pi]), the vertical 
levitation force is determined by the radial magnetic field 
for the PM and the current density of the HTS. The ra- 
dial magnetic field is strongest close to x/a — ±i?pM at 
the bottom of the disk [see Eq. (H)]. Therefore, the ra- 
dial magnetic field and the current close to x/a = ±i?PM 
on the bottom of the disk contribute more to the levi- 
tation force. When the applied field is decreased by a 
small amount, the current close to x/a = ±i?pM on the 
bottom of the disk is reversed to positive first, although 
the total current is still negative (the magnetic moment 
is negative), the levitation force decreases sharply to a 
negative value, because the attractive force results from 
the current close to x/a = ±i?pM at the bottom of the 
disk is much larger than the repulsive force resulting from 
the negative current in the HTS. For a small tr, because 




of relaxation, the current is small, and only a small de- 
crease in the applied magnetic field may lead to a com- 
plete reversal of the current, see Fig. ^, therefore the 
maximum attractive force is closer to zero PM-HTS sep- 
aration. When <T is increased, the current is larger, more 
decrease of the applied magnetic field is needed to reverse 
the current. Therefore the maximum attractive force is 
larger and shifts to larger PM-HTS separation. Further 
increasing a shifts the maximum attractive force to larger 
PM-HTS separation, but the maximum attractive force 
is smaller. This is because much more decrease of the ap- 
plied magnetic field is needed to reverse the larger current 
when a is increased, although the current is larger but 
the radial field is smaller, resulting in a smaller maximum 
attractive force observed in Fig. |6p. 

Another characteristics of the HTS which drastically 
influences the magnetic levitation force is the critical cur- 
rent density of the HTS. The calculated results of ver- 
sus z curves at 6/a = 1.4, cr = 20 and w = 0.1 for different 
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FIG. 5: (a) The vertical magnetic levitation force versus 
the distance 2;atcr = 2,Lj = 0.1 and Jc — 1 for different thick- 
nesses of the superconducting disk. Inset shows the maximum 
repulsive force as a function of the thickness, solid lines are 
guide for the eyes only, (b) The corresponding magnetization 
hysteresis loops. 



FIG. 6: (a) The vertical magnetic levitation force versus 
the distance z at b/a = 1.4, uj — 0.1 and Jc — 1 for dif- 
ferent pinning potentials of the superconducting disk. Inset 
shows the maximum repulsive force as a function of a, the 
solid line is a guide for the eyes only, (b) The corresponding 
magnetization hysteresis loops. 



critical current densities Jc are plotted in Fig. 0a. The 
corresponding magnetization hysteresis loops are shown 
in Fig. For small critical current densities (Jc < 10), 
the calculated results are typical versus z curves, 
which show both the repulsive force branch and the at- 
tractive force branch. The attractive force branches of 
the curves are hidden by the large vertical axis scale. 
However, when Jc is larger than 20, no attractive force 
is observed. When Jc is larger than 60, the F^ versus z 
curves are almost reversible (cr = 80,100). This result 
can also be seen from the magnetization hysteresis loops, 
as shown in Fig. ^a. Increasing the critical current den- 
sity leads to larger hysteresis loop and then to reversible 
magnetization hysteresis loops. Experimentally, F^, ver- 
sus z curves with very srnall hysteresis width has been 
observed in a MTG samplai3, indicating a very high crit- 
ical current density in this sample. 

Based on the critical state model, the constant-field- 



gradient modelEj predicted that the maximum repulsive 
force depends linearly on the critical current density. In 
order to study this dependence, we plot in the inset of 
Fig. ^ the dependence of the maximum repulsive force 
as a function of the critical current density shown as open 
circles. It is obvious that the dependence is a nonlinear 
function. The linear dependence is observed only at very 
low critical current density Jc < 5. A fitting to the ob- 
tained data results in the dependence 

9.3 + Jc ^ ' 

shown as solid line in the inset. From Eq. ( |l^ ) we can 
see that the maximum repulsive force saturates at high 
critical current density. 

Experimentally it is very difficult to get systematic re- 
sults as shown in Fig. ^ and Fig. however, Fig. |6| and 
Fig. M can be used as a quick reference for testing the 
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FIG. 7: (a)The vertical magnetic levitation force versus 
the distance z at b/a = 1.4, ct = 20 and uj = 0.1 for dif- 
ferent critical current densities of the superconducting disk. 
Inset shows the maximum repulsive force as a function of Jc, 
the solid line is a fit with = 3.31Jc/(9.3 + Jc). (b)The 
corresponding magnetization hysteresis loops. 



properties of the sample under investigation. By com- 
paring the experimental results with Fig. ^ and Fig. 
one may have an idea about the pinning potential and 
critical current of the sample. 



D. Force creep 

It has been well known that the current density in HTS 
experiences relaxation (decay with time) , which has been 
extensively studied both theoretically and experimentally 
in terms of the magnetization. Because the vertical mag- 
netic levitation force is determined by the current density 
in the HTS and radial magnetic field generated by the 
PM, it is expected that the levitation force may exhibit 
relaxation, resulting in a time dependent levitation force 
in experiments. However, reports on this dynamic behav- 
ior are very fewBEj. On the other hand, no models so 



far have taken this relaxation behavior of the magnetic 
levitation force into account. This is quite surprising 
considering the operating temperature of the levitation 
system between a HTS and a PM, which is normally at 
the relatively high temperature of 77K, where the relax- 
ation rate is usually quite high. It is also surprising if one 
considers the effects of the force decay on the design of 
any practical applications invoking PM-HTS levitation. 

In this calculation, the force relaxation is taken into ac- 
count by using a voltage-current law E{J) = Ec{J / Jc)^ ■ 
With 1 < n < oo, this material law describes the flux 
creep in terms of the magnetization inside the supercon- 
ducting disk. In Fig. ^ we show the F^. versus z curves at 
b/a — 1.4, a — 2 and Jc = 1 for different frequencies lo at 
which the PM approaches and recedes from the HTS. It 
is obvious that the levitation force depends on the speed 
at which the PM approaches and recedes from the HTS; 
larger speed leads to a larger force loop. On the contrary, 
previous models based on the critical state model assume 
that the levitation force is independent of the speed_at 
which the PM approaches and recedes from the HTSEj. 

In the inset of Fig. |^a, we show the maximum repulsive 
force as a function of the time t = l/uj as open circles 
in a double logarithmic plot. The linear dependence is 
obvious. Fitting to the data gives a, F^ (x t~"^ depen- 
dence of Fz on time. This result actually arises from the 
thermally activated flux motion in the HTS, because the 
voltage current law E{J) = Ec{J/Jc)" is actually a re- 
sult of the current dependence of the activation energy 
U{J) — Uo log(Jc/J); this can be shown by using the Ar- 
rhenius law E = Bv = Bvq exp[-U{J)/kBT]. With this 
activation energy, the relaxation of the current density 
can be derived as J cx leading to the force relax- 

ation shown in the inset of Fig. [^a. Experimentally, the 
force creep has been observed to be logarithmic in timecJ 
within a narrow time window. It is expected that a non- 
logarithmic force creep will be observed if an extended 
time window is measured. 

The force relaxation can also be clearly seen from 
Fig. Ijb, where the magnetization hysteresis loops cor- 
responding to Fig. I^a are plotted. The magnetization 
hysteresis loops increase with increasing frequency. This 
behavior is similar to the so-called dynamic magnetic mo- 
ment relaxation measurement, in which the DC magne- 
tization hysteresis loops are measured at different sweep 
rates of the applied magnetic moment. The magnetiza- 
tion hysteresis loop increases with increasing sweep rate. 
It has been proven that the dynamic magnetic moment 
relaxation is equivalent to the normal relaxation mea- 
surement, in which the magnetic moment is recorded as 
a function of time at fixed temperature aoiLpagnetic 
field, when studying the relaxation of HTScjO. There- 
fore, by measuring the vertical levitation force F^ versus 
z at different frequencies with which the PM approaches 
and recedes away from the HTS, one can study the re- 
laxation of the force as well as the current density of the 
HTS. 
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FIG. 8: (a)The vertical magnetic levitation force versus 
the distance z at h/a = 1.4, o = 2 and Jc = 1 for different 
frequencies at which the PM approaches and recedes from 
the HTS. Inset shows the relaxation of the levitation force. 
(b)The corresponding magnetization hysteresis loops. 



E. Minor force loops and stiffness 

One of the most important parameters used to charac- 
terize the magnetic levitation system using a PM and a 
HTS, is the magnetic stiffness defined as — —dF^/dz, 
which represents the spring constant associated with vi- 
brational motion of a levitation system. The magnetic 
stiffness can be determined by measuring or computing 
minor force loops. The calculated minor force loops at 
different PM-HTS separations for b/a = 1.4, a = 100, 
w = 0.1 and Jc — \ are shown in Fig. ^. Here we use 
an amplitude of 8z = 0.02a, resulting in reversible mi- 
nor force loops. Higher amplitude will lead to hysteretic 
behavior not shown here. The calculated vertical stiff- 
ness is plotted in the inset as a function of the PM-HTS 
distance z. It can be seen from Fig. || that i^T^ > 0, indi- 
cating a stable levitation between a PM and a HTS. 
decreases with increasing PM-HTS separation z. may 




FIG. 9: (a)The vertical magnetic levitation force versus 
the distance z at b/a = 1.4, a — 100, uj — 0.1 and Jc = 1. 
The minor force loops at different distances z are shown with a 
amplitude Sz = 0.02a. Inset shows the stiffness of the system. 



also depend on other parameters, such as Jc, uj, a and 
the shapes of both the PM and the HTS. All these cases 
are easily computed by our numerical method. Detailed 
results will be presented elsewhere. 



IV. CONCLUSIONS 

The current density and magnetic field profiles of a su- 
perconducting disk (SC) with radius a and thickness 2b 
immersed in the non-uniform magnetic field generated by 
a permanent magnet (PM) are calculated from first prin- 
ciples for the superconductor. From the derived current 
density, the magnetic levitation force between the SC and 
the PM has been derived by assuming a voltage-current 
law E{J) = Ec{J/ Jc)^ and a material law B = /ioH. 
The geometry and characteristics of the SC drastically 
influence the vertical levitation force F^. F^ depends 
non-linearly on the critical current density jc of the SC: 
Fz = 3. 31Jc/(9. 3-1- Jc). For thin samples, F^. depends lin- 
early on the thickness of the sample, but beyond a certain 
thickness, F^. is mostly independent of the thickness. The 
flux creep also plays an important role on the magnetic 
levitation force, which is observed to be non-logarithmic 
in time. The stiffness for the SC and PM system has also 
been derived. The calculated results reproduce many of 
the features for the magnetic levitation between a SC and 
a PM exhibited by experiments. 
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